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The theory of optical frequency doubling conversion efficiency is analyzed for the 
small signal input case along with the strong signal depleted input case . Angle phase 
matching and beam focus spot size are discussed and design trades are described which 
maximize conversion efficiency. Experimental conversion efficiencies from the literature , 
which are less than theoretical results at higher input intensities due to saturation , recon- 
version, and higher-order processes, are applied to a case study of an optical communica- 
tions link from Saturn . Double pass conversion efficiencies as high as 45 percent are 
expected. It is believed that even higher conversion efficiencies can be obtained using 
multipass conversion. 


I. Introduction 

In a deep space optical communications link, a frequency 
doubled Nd. YAG laser transmitter emitting at 0.532 pm is the 
preferred transmitter. When the Nd:YAG laser is diode pumped 
the laser is completely solid state and has higher reliability for 
spacecraft communication than its gas and dye counterparts. 
The Nd:YAG laser fundamentally outputs at 1.06 pm. Since 
detectors typically have very low quantum efficiencies at 
1.06 pm, the laser is frequency doubled to 0.532 pm to pro- 
duce a wavelength that has higher detector quantum efficien- 
cies. In this article, frequency doubling conversion efficiencies 
are discussed. First, in Section II, theoretical conversion effi- 
ciencies are discussed. Then, in Section III, experimental 
results from the literature are compared to the theoretical 
results presented in Section II. A case study is then presented 
that uses experimental data from the literature and typical 


optical communication link parameters for a link from Saturn. 
Finally, the conclusion is given in Section IV. 

II. Theoretical Analysis 

In 1961, second harmonic generation was first experimen- 
tally demonstrated with a conversion efficiency of approxi- 
mately 10~ 8 [1]. Today conversion efficiencies of 30 to 
40 percent are not uncommon. Second harmonic generation 
conversion efficiency for plane waves can be expressed as: 
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for small input intensities (see Appendix A [2]). Here, P(2u>)/ 
P(oj) is the conversion efficiency, r?, d Qf{ is the effective non- 
linear coefficient and is the fundamental frequency. The 
term /^ou) denotes the input power while P( 2oo) denotes the 
second harmonic power. The variable L is the crystal length 
and A k equals the difference between the wave vector (k= 2 tt/X), 
k 2 , for the second harmonic frequency, 2co, and the wave 
vector k 1 for the fundamental frequency, co. The beam waist 
radius is co 0 and n is the index of refraction within the crystal. 
The quantities e Q and ju 0 are the electric and magnetic permea- 
bilities, respectively, within a vacuum. 

In Eq. (1), the sin 2 ( AkL/2)/(AkL/2 ) 2 factor can be made 
to approach unity, thus maximizing the conversion efficiency, 
with the proper phase matching techniques. Therefore, an 
understanding of phase matching is in order. If a second har- 
monic wave is produced at one plane and propagates to a sec- 
ond plane, and if another second harmonic wave is generated 
independently at this second plane, then the two waves inter- 
fere. The sin 2 (AkL/2)/(AkL/2) 2 factor describes that inter- 
ference, where Ak = k 2 - 2k x . For maximum efficient fre- 
quency doubling, Ak = 0. Therefore, k 2 ~2k x .A more detailed 
discussion of phase matching is given in Appendix B. 

Looking at our initial conversion efficiency formula, Eq. (1), 
we have an inverse relationship between conversion efficiency 
and beam cross sectional area, go 2 7t, assuming a uniform beam. 
This might lead one to tightly focus the beam. But because a 
tight focus results in greater divergence away from the beam 
waist, as Fig. 1 shows, minimizing the beam waist does not 
necessarily maximize the conversion efficiency. The confocal 
parameter, z, is defined as the distance from the plane of the 
beam waist to the plane in which the beam’s cross sectional 
area is twice that at the beam waist (see Fig. 1). 
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It should be noted that some authors define the confocal 
parameter as twice that distance [3] . 
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This formulation exhibits a conversion efficiency proportional 
to the crystal length, whereas in the plane wave model the con- 
version efficiency is proportional to L 2 . 

An exact analysis by Boyd and Kleinman [4] found that 
the conversion efficiency versus beam -waist relationship fol- 
lows the curve shown in Fig. 2. The optimum conversion re- 
sults when L = 5.68 z. This corresponds to an efficiency 20 per- 
cent greater than if L = 2z. Boyd and Kleinman’s analysis 
incorporated an efficiency reduction factor in the conversion 
efficiency formula to analytically include the effect of focus- 
ing within the crystal. The efficiency reduction factor, h , is a 
function of crystal length, L , and the beam’s confocal param- 
eter, z. When h(L/2z ) enters the small input formula, Eq. (1), 
the yield is [4] : 

In Eq. 4 we assume complete phase matching (sine 2 (A&Z/2) = 1 ) 
and observe that the conversion efficiency is proportional to 
h(L/2z). Therefore, conversion efficiency is maximized with 
the maximum h. Figure 2 illustrates the h(Lj2z ) function ver- 
sus L/2z. Here, and in Eq. (4), we have assumed the most effi- 
cient case where Boyd and Kleinman’s double refraction 
parameter, which includes the effect of the walk-off angle 
between the Poynting vectors of the fundamental and the sec- 
ond harmonic waves [3] , has been set to zero. 





If the confocal parameter is much larger than the crystal 
length, then the beam’s cross sectional area is relatively con- 
stant through the crystal. The plane wave formulation approxi- 
mates this case. However, if the crystal length is much greater 
than 2z, the decreased intensity away from the beam waist 
results in a lower overall conversion efficiency and a focused 
wave formulation must be used. 

If L=2z, we have confocal focusing. By using L=2z=2nncollX, 
area = 7 ro; 2 , and X = 2ttc/lo we obtain an equation for frequency 
doubling conversion efficiency for the focused case : 


All the previous formulations required that conversion to 
the second harmonic be small. In other words, the amount of 
fundamental light traversing the crystal available for conver- 
sion to the second harmonic remained constant throughout 
the crystal. Since an optical communications application re- 
quires higher power levels with much greater conversion effi- 
ciencies, a discussion of the depleted input situation is war- 
ranted. In the depleted input case conversion efficiency is 
high; therefore, the amount of fundamental light traversing 
the crystal is continually reduced; hence, the term depleted 
input. 
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The equation for conversion efficiency for the depleted 
input case is given in Eq. (5) (see Appendix C): 


n 2w) _ 
P( W) 


2 tanh 2 
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Single pass conversion efficiencies for the small signal input 
and depleted input formula are plotted in Fig. 3 for compari- 
son. In Fig. 3, the upper curve is the small signal input formu- 
lation and the lower curve is the depleted input formulation 
for frequency doubling of 1064 nm radiation in a 5.1 mm KTP 
crystal. As can be seen from Fig. 3 the small signal formula is 
only valid for a small region, up to a conversion efficiency of 
approximately 15 percent. Beyond that region, the loss of 
input intensity due to conversion to the second harmonic is 
substantial and the small signal formula breaks down. The 
intensity of the beam decreases as the beam propagates through 
the crystal; thus, the conversion efficiency in the crystal 
decreases. Equation (5) assumes plane wave propagation within 
the crystal; however, the analysis by Boyd and Kleinman for 
focusing within the crystal could be applied to the depleted 
input case as well. Since an optical communication application 
requires high conversion efficiencies, the depleted input situa- 
tion is appropriate. 


TEM 00 operation. Due to the difference in theory versus 
experiment, we use Driscolfs experimental data to calculate 
conversion efficiency in the following case study. 

A specific case for frequency doubling is a cavity dumping 
scheme intended for an M-ary pulse position modulated (PPM) 
optical communications link [6]. As an example, assume a 
pulse repetition rate of 14.3 kHz (114 kbits/s for M=256) with 
a dead time of 44.4 /is and a pulse width of 20 ns. These 
parameters are typical of links being considered from Saturn 
to Earth. If an average laser power of 500 mW is assumed then 
the peak power is 1.75 kW with 35 /ijoules as the correspond- 
ing pulse energy. An average spot diameter of 100 /im inside 
the crystal corresponds to a power density in the crystal of 
22 MW/cm 2 . Driscolfs graph of experimental data [5] for 
double pass efficiency through a 5.1 mm KTP crystal relates a 
22 MW/cm 2 beam to approximately 28 percent double pass 
efficiency as shown in Fig. 5 with an arrow. 

That conversion efficiency can be improved upon by de- 
creasing the repetition rate, decreasing the average spot size, 
or increasing the average power level. Table 1 summarizes the 
results for various system parameters. Note that 17 in Table 1 
only results from a double pass senario [5] and will probably 
be higher for multipass configurations. 


III. Case Study 

Let us examine a specific nonlinear frequency doubling 
crystal, KTiOP0 4 (KTP). The KTP crystal is a relatively new 
second harmonic generation material and has a few advantages 
over such predecessors as KDP and LiNb0 3 . It has a high dam- 
age threshold and a high nonlinear coefficient, as well as 
excellent optical quality. 

The d eff values within the literature vary between theoreti- 
cal and experimental results. By averaging the values of <2 eff , 
which were determined from the graphical data presented in a 
Lockheed study of KTP [5] , a value of 2.2 X 10“ 23 m/V was 
found. In Fig. 4 experimental values for single pass second- 
harmonic-generation conversion efficiency for a guassian beam 
have been graphed with theoretical values calculated in Sec- 
tion II for the depleted input case. At fundamental intensities 
less than 80 MW/cm 2 the two curves are in agreement. How- 
ever, for intensities greater than 80 MW/cm 2 the experimental 
values are less than the expected theoretical results. Driscoll 
et al attribute this saturation of reconversion to the fundamen- 
tal and higher-order processes [5]. Driscoll et al achieved 
higher conversion efficiencies with a multimode laser, proba- 
bly due to localized higher intensity within the modal spot 
structure; however, for our application of deep space optical 
communication, we are primarily interested in single mode 


A better understanding of the roles of spot size and pulse 
width can be gained by breaking down Table 1 and graphing 
the information. Figure 6 plots average power in watts versus 
double pass conversion efficiency for different spot sizes. Fig- 
ure 6 shows that for the 50 jum beam, the conversion effi- 
ciency peaks and begins to drop off at approximately 500 mW, 
while the beam with a spot size of 100 /im peaks and begins to 
drop off at 2 watts. For two beams of the same average power, 
halving the spot diameter quadruples the intensity. At high 
intensities, the conversion efficiency decreases due to satura- 
tion, reconversion, higher order processes, and intensity fluc- 
tuations [5] . Thus, the curve for the 50 jum beam peaks and 
begins to dip at average power levels which are one fourth of 
those for similar efficiencies with a 100 /im beam. Figure 7 dis- 
plays curves of double pass conversion efficiency versus aver- 
age input power with a 10 ns pulse width and a 20 ns pulse 
width. For the specific case calculated, conversion efficiency 
peaks at 1.25 watts for the 10 ns pulse width case at a nominal 
pulse repetition rate of 14.3 kHz as discussed earlier. The 
20 ns pulse width case peaks at a slightly higher average power 
level of 2.0 watts for the same repetition rate. 

Furthermore, by varying the parameters, as the graphs and 
table display, double pass conversion efficiencies from 28 per- 
cent to 45 percent should be attainable for an optical com- 
munications link. However, even higher efficiencies are believed 
possible with a multipass system. 
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levels in this case, for a uniform average 100 pm diameter spot 
size, were approximately 90 MW/cm 2 which was well below 
the damage threshold of ~350 MW/cm 2 . Experimental conver- 
sion efficiencies found in the literature at this intensity level 
were in close agreement with those calculated from the theo- 
retical discussion presented in Section I. At higher intensity 
levels, experimental results were less than those calculated due 
to saturation from reconversion and higher-order processes. 
With a multipass configuration, it is possible that still higher 
conversion efficiencies may be achievable. 
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IV. Conclusion 

A theoretical analysis of frequency doubling conversion 
efficiency has been presented along with a case study of an 
optical communications link from Saturn. The case study used 
experimental test results found in the literature and applied 
them to typical link parameters. An average power, 1 watt, 
laser operating at over 100 kbits/s in a pulse position modula- 
tion mode with M=256, and a pulse width of 10 ns, yielded 
double pass conversion efficiencies as high as 45 percent. Flux 
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Table 1. Summary of various system parameters 


Spot 

diameter 

Pulse 

width 

Dead 

time 

Rep rate 
M=256 

Average 

power 

Peak 

power 

Power 

density 

n 

100 Mm 

20 ns 

44.4 ms 

14.3 kHz 

500 mW 

1.8 kW 

22 MW/cm 2 

28% 

100 Mm 

20 ns 

44.4 ms 

14.3 kHz 

1 W 

3.5 kW 

45 MW/cm 2 

40% 

100 Mm 

20 ns 

44.4 ms 

14.3 kHz 

2 W 

7.0 kW 

89 MW/cm 2 

45% 

100 Mm 

20 ns 

44.4 ms 

14.3 kHz 

3 W 

10.5 kW 

134 MW/cm 2 

43% 

50 Mm 

20 ns 

44.4 ms 

14.3 kHz 

250 mW 

0.88 kW 

45 MW/cm 2 

41% 

50 Mm 

20 ns 

44.4 ms 

14.3 kHz 

500 mW 

1.8 kW 

89 MW/cm 2 

45% 

50 Mm 

20 ns 

44.4 ms 

14.3 kHz 

1 W 

3.5 kW 

179 MW/cm 2 

43% 

50 Mm 

20 ns 

44.4 ms 

14.3 kHz 

2 W 

7.0 kW 

357 MW/cm 2 * 

40% 

100 Mm 

10 ns 

44.4 ms 

14.3 kHz 

500 mW 

3.5 kW 

45 MW/cm 2 

40% 

100 Mm 

10 ns 

44.4 ms 

14.3 kHz 

1 W 

7.0 kW 

89 MW/cm 2 

45% 

100 Mm 

10 ns 

44.4 ms 

14.3 kHz 

2 W 

14.0 kW 

178 MW/cm 2 

43% 

^damaged threshold 







FUNDAMENTAL INTENSITY, MW/cm 2 


Fig. 4. Theoretical and experimental frequency doubling conver- 
sion efficiencies are plotted for single pass through a 5.1 mm KTP 
crystal. The saturation of the experimental results is due to recon- 
version and higher order processes. 



AVERAGE POWER, watts 


Fig. 6. Conversion efficiency versus average laser power for a PPM 
optical communications link. Assumed data rate is ~-100 kbits with 
a pulse width of 20 ns and a dead time of 44.4 fis. Efficiencies for a 50 
M m and 100 ^m average spot size are graphed. 
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FUNDAMENTAL INTENSITY, MW/cm 2 


Fig. 5. Experimental conversion efficiency versus incident power 
density for a Gaussian 1.06 fj.m fundamental beam. The crystal was 
5.1 mm in length [5]. 



AVERAGE POWER, watts 

Fig. 7. Conversion efficiency versus average laser power for a PPM 
optical communications link. Assumed data rate is — 1 00 kbits with 
44.4 (jlS dead time. Average spot size in the crystal was 100 ^m. 
Efficiencies for pulse widths of 10 ns and 20 ns are graphed in the 
figure. 
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Appendix A 
Conversion Efficiency 


Crystals capable of higher order processes have nonlinear 
terms included in their polarization. 

P = bE{ 1 + a x E + a 2 E i 2 + a 3 E 3 + . . .) 

where a { and b are constants and E is the electric field. Speci- 
fically, let us examine the effect of the first nonlinearity, 
P = dE 2 y which is responsible for second harmonic generation. 

Starting with a one dimensional electric field traveling in 
the z direction, let us consider traveling waves of three fre- 
quencies, coj , co 2 , and co 3 . 


k ' 2 ’ ’ ' 2 1 2k 

E.^,2,1) =I[£ 3 .( 

Then, using the equation 


a2 e = a It + M ° e 


e i^ l t-k 1 z) + cc] 

(Al) 

)e 'w-M) +ca] 

(A2) 

e ^3 f - fc 3 Z ) + c . c .J 

(A3) 

9£ 2 92 

+ n 

(A4) 

b t 2 0 




d eff E U E lk 


e -i{k x +k 2 -k 2 )z 


(A7) 


For second harmonic generation, co l = co 2 and co 3 = co x + 
CO 2 = 2 OJj . 


If we assume that the quantity of power lost by the input 
beam at co i due to second harmonic conversion is negligible, 
then 


As a result of the above two conditions only Eq. (A7) need be 
examined. Furthermore, the medium chosen will usually be 
transparent to radiation at c o 3 , hence a 3 = 0. 


Therefore, Eq. (A7) reduces to the following: 




^eff E \ i k 


iAkz 

e 


(A8) 


where 


where (P NL ) t = d eff EjE k is the nonlinear polarization and a 
is the conductivity, it can be shown that 



i co. 


^ ‘ d eff E 3/ E 2k + e 




(A5) 


Ak = k 3 ~k 1 ~k 1 

If we take E 3/ (0) = 0 for no second harmonic input and let 
the. frequency doubler’s length be L, then the solution to Eq. 
(A8) is 

nr e tAkL _ j 

■ -“/t d „, E u E n — SF“ (A9> 

where multiplication by E 3/ . leads to 





E U E 3/ 


e -i(k 1 -k 3 + k 2 )z 


(A6) 


e 3 .(de;.(l ) = 4 w 2 d 


eff 


x e 2 .e 2 l 2 

1 1 Ik 



(A 10) 
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This equation leads to an expression in terms of second har- Utilizing a similar expression for / > (co)/Area, an expression 

monic output power when for conversion efficiency, P(2u))/P(u>) can be found to be: 


?( 2<^) _ i r±_ F p* 

Area 2 J u 3 / c 3j 


(All) 


is substituted into it. Thus: 

P(2co) _ 1 JK 2 (J \2 p2 p2 r 2 
“Area - " 2x/T W E U E Vc L 


sin 


AkL 




(A 12) 


n 2a,) jk\ , n (tt) 

” ' " Vo) A '« fAkL \ ! 

(A 13) 

(note that e = e 3 and e x = e 3 = e Q rt 2 ). See main text for 
definition of variables. 
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Appendix B 

Crystal Angle Phase Matching 


If the crystal is birefringent, having axially dependent indi- 
ces of refraction, then k(<uo) = c oQie^ 1 ^ 2 n( to). The k condi- 
tion of k 2 = 2 k x , as discussed in the main text, translates to an 
«(2co) = «(co) condition. Since in dispersive materials n is pro- 
portional to co, it is not possible to match n{2co) and n( co) if 
the wave is propagating down a single crystal axis, because 
h(2co) will always be larger. Therefore, in order to satisfy the 
equality of /?(co) = a?(2co), a combination of the crystal indi- 
ces of refraction is necessary. For uniaxial crystals, two axis 
indices are the same: n x = n y - n Q and n z = n e where z is the 
crystal optic axis, n Q is the ordinary index of refraction and n e 
is the extraordinary index of refraction. Figure B1 is a graphi- 
cal representation summarizing the indices of refraction for 
two orthogonal polarizations at varying angular propagation 
directions to the optic axis for differing wavelengths. In the 
graphs we are considering the case of the negative uniaxial 
crystal, n Q > n e . The wave-vector, D 0 , propagates along s with 
an electric field polarized normal to the page. Thus, D 0 
remains constant, as seen by the circular representation of 
for varying propagation angles, 6 . On the other hand, D e has 


an electric field polarized perpendicular to D 0 and the direc- 
tion of propagation; thus, its magnitude varies for different 
propagation angles, 6. Figure Bl(a) illustrates these index rela- 
tionships at frequency co while Fig. Bl(b) illustrates them at 
2co. Note the increased size at 2to due to the higher frequency. 
In order to obtain n(co) - «(2co), the 2 curves are overlapped. 
The intersection of the two curves yields the desired answer. 
The solution results from using n Q at oo and n e (Q) at 2co as 
seen in Fig. Bl(c). This angular solution is analytically described 
as follows: 

1 _ cos 2 0 sin 2 # . 

' ■ ■ ^ / 

(n e (2<jj,d)) 2 (n 0 (oj)) 2 (n e (2oo)) 2 

where 6 is the angle between the ray’s axis of propagation and 
the optic axis. This same type of analysis is used with biaxial 
crystals where n x ¥ i n y ^n z . Since all axes have different 
indices of refraction, the analysis is much more complicated. 
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Fig. B1. Phase matching curves for ordinary and extraordinary indices of refraction 




Appendix C 

Depleted Input Conversion Efficiency 


Let us utilize Eqs. (A5), (A6), and (A7) from Appendix A 
and follow a similar approach as found in Yariv [2] . 


dE u = 

. - 0 ^ 

£*, . 

-l (, 0 ) 1/2 

dz 


1 * 

2(e 1 ) 1 / 2 

dE l k _ 


#r* 

+ - 2 (, 0 ) 1/2 

dz 

2(e 2 y' 2 


2(e 2 y/ 2 

dE v . 

- CT 3 (M 0 y /2 

p 

- 3 m /2 

dz 

2(e 3 ) 1/2 

*3 i 

2(e 3 ) 1/2 


d .« E 3i Kk*- iAkz 


d „E, .E*.e iAkz 

eff 1 1 3/ 


d „£..£* e iAkz 
eff 1/ 2k 


We can set o t = o 2 = o 3 = 0 because the medium chosen 
should be transparent to radiation of the input frequency and 
the second harmonic. Furthermore, for frequency doubling, 
E x = E 2 . Since the phases should be matched, Ak = 0. Also, 
if £^(0) is real, then E x (z) is real. Let us drop the ///c notation 
and define E 3 = - iE 3 . 

Therefore, Eqs. (Cl) and (C3) become: 
dE —co (ii yt 2 

1 1 0 t r?> T7* try a\ 




d K 1/2 




:f— 2 — rf ff £? 

2 ( y/2 err , 


Therefore, 


— £ 2 + E'} 


, 2 w i( e x ) 1/2 


w,(ej 1/2 


If there is no input at co 3 , then integrating the above equa- 
tion and substituting e Q n 2 = e r « e 3 yields: 


(co ) 

E ? +E ? (4) = ^ (0) 


Substituting into Eq. (C5) leads to: 
co>J 1/2 / 

3 _ 3 O , I rr 2 tr\\ 




(g « °) S 2 ) 


Integration of Eq. (C8) leads to: 

_ / w a\ 1/2 . , Kff^o^i o J3 ) 1/2 £i(0). 


£3 = *,«» 


2 (e G « 3 ) 1/2 


Given that 

I. P 

where A - beam area, and using co 3 = 2cOj , Eq. (C9) becomes: 




|£, ( 0 ) | 2 


^l/ 2 £ 3/4 „3/2 


Making the small angle approximation leads to the small signal 
formula of Appendix A for conversion efficiency. 
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